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Abstract 

Building on mathematical similarities between quantum mechan¬ 
ics and theories of diffusion-influenced reactions, we discuss how the 
propagator of a reacting molecule pair can be represented as a product 
of three factors in the Laplace domain. This representation offers sev¬ 
eral advantages. First, the full propagator can be calculated without 
ever having to solve the corresponding partial differential equation or 
path integral. Second, the representation is quite general and capa¬ 
ble of capturing not only the classical Smoluchowski-Collins-Kimball 
model, but also alternative theories, as is here exemplified by the case 
of a delta- and step-function potential in one and two dimensions, re¬ 
spectively. Third, the three factors correspond to physical quantities 
that feature prominently in stochastic spatially-resolved simulation al¬ 
gorithms and hence the interpretation of current and the design of 
future algorithms may benefit. Finally, the representation may serve 
as a suitable starting point for numerical approximations that could 
be employed to enhance the efficiency of stochastic simulations. 


1 Introduction 

Stochastic processes, like Brownian motion, and quantum phenomena can 
be described by mathematical structures that appear to be quite similar 
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MUM- A case in point is the transition probability density function 
(PDF) p(x,t\xo,to), a central quantity in the theory of Brownian motion 
that yields the probability to find a Brownian walker at x at time f, provided 
it started at xq at to < t. The transition PDF can be calculated in terms of 
a path integral that may formally be written as 

p(x,t\x 0 ,t 0 ) ~ ^2 e -Ss[s(T)] ; (1) 

x(t):x o~>x 

where the sum includes all paths x(t) that satisfy x(to) = xo and x(t) = 
x. 5^[a;(r)] denotes a certain functional of x(t). Quite analogously, the 
Feynmann transition amplitude K(x, t\xo, to) is also given by a path integral 

K(x,t\x 0 ,t 0 ) ~ ^2 eH s[x(r)1 , (2) 

x(t):x q—>x 


where <S[x(t)] refers to the classical action functional. Although the physical 
context decribed by Eq. m is quite different from the one represented by 
Eq. ([2]) . the path integrals can in general be transformed into each other by 
means of a Wick rotation [2T1 l6l [22] and the functional Sb[x(t)] in Eq. (P) 
is the Euclidean version of the action <5[x(r)]. 

Typically, path integrals come along with partial differential equations 
(PDE) and integral equations [26] . 

Hence, unsurprisingly, the formal correspondence continues on the level 
of the associated PDE. As is well-known, both the transition PDF p(x, t\xo, to) 
and the Feynmann amplitude K(x,t\xo,to) are the Green’s functions (GF) 
of the diffusion equation 


dp(x,t\x 0 ) 

dt 

P(x,t = 0|x 0 ) 


8 2 

6(x - x 0 ), 


p(x,t\x o), 


( 3 ) 

( 4 ) 


and the Schrodinger equation, 


dK(x , t\xo) 

Ft 

K(x,t = 0|xq) 


' in d 2 

2m dx 2 
5(x - x 0 ), 



K(x,t\x 0 ), 


( 5 ) 

( 6 ) 


respectively. 

In this manuscript, we exploit the fundamental relationship between 
quantum mechanics (QM) and Brownian motion in the context of theories of 
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diffusion-influenced reactions [12, 20]. In the classical Smoluchowski-Collins- 
Kimball (SCK) picture [251171 [T2l 120] . the diffusion-influenced reaction of an 
isolated pair of molecules is decribed by a diffusion equation (Eq. (|3|l) with 
vanishing potential V ( x ). The actual reaction is instead implemented by im¬ 
posing suitable boundary conditions (BC) at the so-called encounter radius 
However, there is an alternative approach that abandons 
the BC and postulates a non-vanishing potential I27isi uni m m nans]- 
This formulation is not only capable of including the classical theory, but 
it may also serve as a starting point to construct more general theories of 
diffusion-influenced reactions. At the same time, its mathematical structure 
is closer to the formalism of QM and hence it may pave the way to borrow 
techniques from QM for theories of diffusion-influenced reactions. 

The manuscript is structured as follows. First, we briefly describe the 
path decomposition expansion (PDX) that has been applied to the quantum 
mechanical propagator. Then, we briefly review the theory of diffusion- 
influenced reactions and discuss the PDX for the classical theory and the 
Doi (also referred to as volume reactivity) model. We explicitly show that 
we can recover the GF of an isolated pair from the PDX for the case of the 
SCK theory (implemented as ^-function potential) in one dimension (ID) 
and for the volume reactivity model in 2D. 


2 Theory 

2.1 The Path Decomposition Expansion 

The PDX was developed in Refs. [1] to study quantum tunneling, subse¬ 
quently it has been used for other applications as well, see Refs. H3E2|. in 
Refs. (2EIE93, this technique was used to calculate the Feynmann propagator 
for a 5- and step-function potential in an elegant and simple way. We point 
out that these potentials correspond exactly to the SCK and Doi model in 
the context of theories of diffusion-influenced reactions, as we will expand 
on later. Here, we consider a quantum mechanical particle in ID. Then, the 
propagator satisfies the following convolution relation [3, 22l [28l f29lj : 


K(x,t\x 0 ) 


iH f 1 T r , . . 

2m Jo K ( X,t \ a,T ^ 


dK r (d,T |so,0) 


di 


dr. 


i=a 


(7) 


This relation means that every path can be split into a restricted part (rep¬ 
resented by K r ) that starts at xq at time to = 0, never crosses the boundary 
at x = a = 0, but that ends at the boundary at time r, and an unrestricted 
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part (represented by K) that starts at the boundary at time r and ends at 
x at time t. Similarly, one has 


K(x,t\x 0 ) 


ih r l dK r (x,t\£,r) 

2m J 0 d£ 


K(a, t\ xo, 0 )dr, 


( 8 ) 


which corresponds to a split into into an unrestricted part that starts at Xo 
at time to = 0, may cross the boundary x = a = 0 many times, and finally 
ends at the boundary after time r and a restricted part that starts at the 
boundary at time r, never crosses the boundary after that and ends at x at 
time t. Combining these two formulas (Eqs. (0, (jH|) ), one obtains for the 
Feynmann propagator 


K(x,t\x 0 ) 


U«L 


4 m 2 
xK(a, T\a, r) 


T dT dK r (x,t\£,T) 


dK r (^,r|x 0 ,0) 


£=a 


5^ 




(9) 


Refs. [28], [29] exploit that the three quantities on the rhs of Eq. ([9|) can 
be easier calculated than the full propagator K(x,t\xo). This convolution 
relation (Eq. Q) is the central identity that we will use in the following and 
we will explicitly check that it holds true in theories of diffusion-influenced 
reactions as well. 


2.2 Smoluchowski-Collins-Kimball theory 

The theory of diffusion-influenced reactions is traditionally described by the 
diffusion equation (Eq. (|3|) ) with vanishing C(x), subject to a suitable BC 
that implements the phyics at the encounter distance. The Collins-Kimball 
BC [7], also referred to as radiation BC, generalizes the Smoluchowski (or 
absorbing) BC |25| and reads 

= KaPrad(x = a, t\xo). (10) 

x—a 

In the limits K a —> oo and K a —> 0, the radiation BC reduces to the absorbing 
and non-reactive (reflective) BC, respectively. The survival probability is 
defined as 

/•OO 

S ra ,d{t) = p Iad (x,t\x 0 )dx. (11) 

J a 


D dp radOMR)) 
dx 
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Now we may invoke the diffusion equation (Eq. (|3|), V(x) = 0) and the 
definition of the survival probability (Eq. fTTl) , which tells us that 


dS iad (t\x 0 ) 

dt 


^ Prad(M|s 0 ) 

dx 


x=a 


( 12 ) 


Note that this relation is true independent of the imposed BC at the en¬ 
counter distance. However, if we impose an absorbing BC, the negative time 
derivative of the survival probability bears a special meaning. It is the first 
passage time PDF 


fpp(t\x 0 ) : = 


5<S' a bs(t|3'o) 

dt 


(13) 


Furthermore, in the case of a radiation BC, the quantity K a Prad(«) t\a) is the 
rebinding time PDF that plays an essential role in understanding spatial 
stochastic fluctuations [Ml HI)] . Owing to the Collins-Kimball BC (Eq. ([TUI) 
and Eq. (fT2l) one has 


frebit) := _ d^ r ad(^° -^ = Ka p rad ( a , £|a). (14) 

at 

Finally, it will turn out to be convenient to consider the PDF /lr(t\x, t) for 
the last reflection time r before t, given that the molecule is located at x at 
time t, which is closely related to the first passage time PDF via 


fLR{r\x,t) = f FP {t - t\x). 


(15) 


As already mentioned, an alternative and more general approach, that 
abandons the requirement of a BC at the encounter distance, implements 
the chemical reaction by adding appropriate sink terms to the equation of 
motion, i.e. one deals with a diffusion equation featuring a non-vanishing 
potential (Eq. (|3|)). For instance, the SCK model may be defined by the 
diffusion equation with V{x) = —K, a 5(x — a), while for volume reaction 
theories the potential assumes the form V(x) = —n r Q{a — x). 

In the following, we will focus on these two examples and show that the 
corresponding full propagator can be obtained by a convolution relation that 
is analog to Eq. ([9j). More precisely, one has 


Pv(x,t |x 0 ) 


±D 


I dT I 

Pv(a, T — r\a) 


T dr dffabsQM ~ T\0 


OH 

dPab S (a,T\x 0 ) 


£=a 


da 




(16) 
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Note that if both x,xo > a or x,x$ < a, one has to add a second term 
—Pabs{x, t|xo) on the lhs of Eq. (fT6]h These two cases also correspond to the 
positive sign on the rhs, while the cases x < a, xq > a and x > a, xo < a 
correspond to the negative sign. Furthermore, Eq. (1161) shows that the 
restricted propagator K r in the quantum mechanical version (Eq. @) cor¬ 
responds to the GF p a bs with absorbing BC. However, we emphasize that 
p a bs may describe propagation in the presence of an absorbing boundary and 
a non-vanishing potential. 

2.3 Partially absorbing trap in one dimension 

To check the validity of Eq. (1161) . we consider a partially absorbing trap in 
one dimension at x = a = 0. Instead of following the classical route via 
imposing a radiation BC, one may employ the following equation of motion 


[HI [23] 



(17) 


The GF solution p$(x,t |xo) has been presented in Ref. [8]. Note that this 
GF is not given by exactly the same expression as the solution p ia d{x, t\xo) 
of the corresponding radiation BC problem, as detailed in Ref. [23]. In 
particular, the GF ps(x,t\xo) is defined for all —oo < x,xo < oo, in contrast 
to the SCK solution p Ta d( x ,t\xo) that is restricted to the semi-infinite line, 
0 < x,xq < oo. In the following, we will focus on the case x < 0, xo > 0, 
but all other cases can be treated similarly. 

We observe that in the (5-function case, the convolution relation (Eq. (1161) 
may be rewritten in the Laplace domain as product of three prominent 
waiting time PDF. To see this, we first note that in the present context one 
has to take into account two different first passage time PDF, corresponding 
to xo < 0 and xq > 0. More precisely, 



Furthermore, one has (see Eqs. (ED, ED, ED) 


fui(T\x,t) = 


fpp(t\xo) 



x<0, (20) 


( 21 ) 
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Next, we would also like to point out that, although p va ,d(x,t\xo) is different 
from ps(x,t |xq) for general x,xq, one has 


Prad(0,t|0) = 2p s (0,t\0), (22) 

and 

Ss(t\x 0 ) = S iad (t\x 0 ), x 0 > 0, (23) 

where 

/•oo /-a=0 

Ss(t\xo):=Sf(t\x Q ) + Sf(t\x Q ):= p 5 (x,t\x 0 )dx+ p 5 (x,t\x 0 )dx 

J a= 0 J— oo 

(24) 

as discussed in Ref. [23]. Thus, using Eqs. (fl6l) . (fTTl) . (fl5l) . f[20]) . (1211) and 
(1221) . we can substitute the first passage and rebinding time PDF for the 
spatial derivative terms and for ps(a, s|a), respectively, in Eq. (fl6l) to obtain 

ps(x,s\x 0 ) = -^—f^ R (s\x)f reh (s)f^p(s\x 0 ), x < 0, xq > 0, (25) 

ZK, a 

where we employed the following notation for the Laplace transform of a 
function g(t ) 

POO 

£[ 9 }{s):=g(s):= e~ st g{t)dt. (26) 

■Jo 

One may also verify Eq. ()25l) explicitly. Transforming the solution given 
in Ref. [8] to the Laplace domain, one arrives at 

e -v(\x\+x 0 ) 

Ps{x,s\x 0 ) = 9 + h y x<0,xo>0, (27) 


where v := ys/D and h := n a /D. Then, using the expressions for the GF 
with absorbing BC in the Laplace domain [5] 


PabsO> S M 

PabsOD S M 

we obtain 


1 

~2Dv 

1 

~2Dv 


e ~v\x~xo\ 

—v\x—xo\ 


e ~v\x+xo\ 

e ~v(x+x 0 ) 


X , Xq < 0 , 
X , Xo > 0 , 


_ ^1^0 

fpp(s\x 0 ) = e V So, 
fpp(s\x 0 ) = e _v ^v^5. 


(28) 

(29) 


(30) 

(31) 


Using Eqs. (l27|) . (l30ll and (13TT) . we see that Eq. (1251) does indeed hold true. 


7 










2.4 Area reactivity model 

The second example that we discuss is the area reactivity model in 2D 
OH UM 119] . also referred to as Doi model in the case of an irreversible 
reaction mm®- The equation of motion reads 


dpe{r,t\r 0 ) 

dt 


D 


8 2 1 _ 3 \ 

dr 2 r dr J 


K r ®(a 


r) 


Pe(r,t\r 0 ). 


(32) 


In the following, we will show that 


P©(Ds|r 0 ) 


±4A 2 D 2a# » b - (r ’ s|{) 


9{ 


p e {a,s\a) 


dpabs(C, ^ko) 


5=« 




C=a 


(33) 


Again, if both r, ro > a or r, ro < a, one has to add a second term — p a bs(L s|ro) 
on the lhs of Eq. (HTH) . as in the ID case discussed before. 

Furthermore, we emphasize that, although we consider here explicitly the 
irreversible reaction, Eq. (1331) remains valid for the reversible case. In fact, 
one only has to substitute w := ^/(s + n r )/s —>• (s + K r + Kd)/(s + , 
see Ref. [T9] . 

To be definite, we consider the case tq > a, r < a. Then, one has mm 


Pe{r,s\r 0 ) 

M 


_ I 0 (wr)K 0 (vr 0 ) 

2iraDJ\f 

= vIo(wa)Ki(va) + wI\(wa)Ko(va), 


(34) 

(35) 


where /o,i, A"o,i refer to the modified Bessel function of first and second kind, 
respectively, and zeroth and first order, [l] Sec. 9.6]. 

Let us now turn to the GF p a bs(L t\ro) satisfying absorbing BC. Similarly 
to the 1-D (5-function case considered before, one has to take into account 
two different GF p a b s , corresponding to the domains r, ro > a and r, ro < a. 
For r > a,ro > a, one has in the Laplace domain 


pI\ bs( r >sM Pfree(r,v\r 0 ) K 0 (vr)K 0 {vr 0 ) \ 

ZTTJJ Kq[VQ) 

(36) 

where 



~ , . ) _ 1 f Io(vr 0 )K 0 (vr), r > r 0 
p {iee [r,s\r 0 ) 27 r D \ Io(vr) K o( Vr o), r < r 0 

(37) 

It follows that 



2 ™-D 8j?b! £ SW 

K 0 (vr 0 ) 

(38) 

dt 

e=a Kq(vo) 
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Next, for r < a,ro < a, one has 


PabsO^M = Phee{r,w\r 0 ) 


1 

2 ^D h 


(■ wr)I 0 (wr 0 ) 


K 0 (wa) 
Io(wa) ' 


(39) 


Note that the appearance of w instead of v is due to the fact that within the 
domain r < a the molecule moves in the presence of a constant potential 
given by — n r . Thus, one obtains 


2iraD 


ftPabsMIO 




Ipjwr) 
Io(wa)' 


(40) 


Finally, using Eqs. (l34l) . (l38l) and (j40]l . we can verify that indeed Eq. (l33l) 
does hold true. Similarly, one can show that Eq. (|33D is fulfilled by the GF 
corresponding to r, > a and r, ro < a, as well as by the GF describing the 
reversible reaction. 

Finally, note that if we used for r < a the expression ~dS< hs (t\r)/dt, 
instead of the flux (Eq. (1401) we would not recover the full propagator. The 
reason is that in the case of a (5-function potential the flux is equal to the neg¬ 
ative time derivative of the survival probability (up to a sign, cf. Eqs. (l20j) . 
(j2T]0 . However, in the AR model, due to the ©-potential, one has instead 

m 


dt 


= 2iraD 


dPabsMM 

dr 


r=a 


KrS< hs (t\r 0 ). 


(41) 
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